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Abstract 

A global analysis of the scaling behaviour of a system with a 
scalar order parameter quenched to zero temperature is obtained 
by numerical simulation of the Ginzburg-Landau equation with 
conserved and non conserved order parameter. A rich structure 
emerges, characterized by early and asymptotic scaling regimes, 
separated by a crossover. The interplay among different dynam- 
ical behaviours is investigated by varying the parameters of the 
quench and can be interpreted as due to the competition of dif- 
ferent dynamical fixed points. 



1 



1 Introduction 



In recent years a great effort has been made to understand the scaling be- 
haviour observed in the late stage of the phase ordering process following 
the temperature quench of a system from an initial disordered state to a 
final state inside the coexistence region, below the critical point. A the- 
oretical framework for the description of this phenomenon, well suited 
both for analytical and numerical approaches, is provided by the time 
dependent Ginzburg-Landau model (TDGL) [[[]. 

In the late stage of phase ordering the order parameter saturates 
locally to the equilibrium values, giving rise to configurations with do- 
mains separated by sharp interfaces. In this regime the only dynamics 
left in the system is, therefore, interface motion and the subsequent time 
evolution is characterized by the coarsening of domains, while their mor- 
phology remains invariant. This leads to the formulation of the scaling 
hypothesis 0, whereby the residual time dependence in the system is 
due only to the growth of the size of domains according to the power 
law L(t) ~ with z = 2 for non conserved order parameter (NCOP) 
and z = 3 in the conserved case (COP). According to this hypothesis, 
the equal time order parameter correlation function obeys an asymptotic 
form of the type 

G{\x-x'\,t) ~L~ a {t)F 
where a = due to the formation of compact domains and F is a scaling 



\x — x 

~1W 



(1) 
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function. This in turn implies the following form 

C(k,t) ~ L{t) d - a F[kL(t)} (2) 

for the structure factor, which is the space Fourier transform of G(\x — 
x'\,t). Although a full theory of scaling is yet to come, the prediction of 
the scaling hypothesis for the correlation function is confirmed by exper- 
iments ||, numerical simulations Q and by exactly soluble models for 
NCOP 0. Less is known about the behaviour of quenched systems in the 
regime preceeding the asymptotic dynamics |J. At early times domains 
are quickly formed but the order parameter is still away from saturation 
inside the ordered regions. As a consequence much more freedom is left 
to the system as compared to the asymptotic regime, since modulations 
of the field are also possible. 

In this paper we present a global analysis of the scaling behaviours 
obtained by simulating the TDGL equation for a two-dimensional system 
and by varying the parameters of the Hamiltonian over the T = mani- 
fold of the equilibrium phase diagram. The time evolution of the system 
is followed from the instant of the quench down to equilibration. In so do- 
ing, we uncover a structure much more rich than usually realized. In the 
very early stage, for t smaller than a crossover time ti, the system relaxes 
towards equilibrium with a purely diffusive behaviour. During this early 
regime the amplitude of the order parameter shrinks to zero in order to 
remove spatial inhomogeneities and the correlation function obeys the 
scaling form (0) with exponents a = d and z = 2 or z = 4 respectively 
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for NCOP or COP. Then, at t ~ t±, the system enters an intermediate 
regime characterized by exponential growth of the order parameter to- 
wards its local equilibrium value. Eventually, after a characteristic time 
tii the late stage scaling is asymptotically obeyed, with exponents a = 
and z = 2 or z = 3 for NCOP or COP. This whole structure is schemat- 
ically represented in fig. 1 with symbols to be specified in the following 
section. We stress that what we call early stage here precedes the usual 
early time behaviour characterized by exponential growth. For common 
choices of the parameters entering the model the crossover time t\ is too 
short to make the early stage observable. However, this can be greatly 
amplified by a proper choice of the parameters of the TDGL equation. 
For COP this amplification can be obtained by performing asymmetric 
(off-critical) quenches. 

In the present paper this structure is investigated by the numerical 
solution of the TDGL equation and is interpreted in terms of the interplay 
between different fixed points of competing stability. For short times 
the dynamics is controlled by the trivial fixed point of simple diffusion 
whereas asymptotically the attractive ordering fixed point dominates. 
The intermediate regime corresponds to the crossover in the time interval 
t\ < t < t2 ■ This interpretation is supported by the comparison with 
the exactly soluble large- N model, where similar properties are observed 
analytically [[7). 

The outline of the paper is as follows: in the next section we intro- 
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duce the TDGL equation and set up the notation. In section 3 we present 
the results of numerical simulations of the model in the scalar case which 
clearly show the existence of the early scaling regime. Data are presented 
for conserved and non-conserved order parameter, critical and off-critical 
quenches. The role of the parameters entering the TDGL equation and 
their effect on the early stage behaviour is elucidated. Section 4 is ded- 
icated to the comparison with the exact solution of the large- iV model. 
Finally, in section 5 a summary of the results is presented and concluding 
remarks are made. 

2 The model 

We consider a system with a scalar order parameter t) initially pre- 
pared in a configuration $(5, 0) sampled from a high temperature uncor- 
rected state with expectations 

<$(£,0)>=0 (3) 

and 

< $(£, 0)$(f, 0) >= AS(x-x') (4) 

where A is a constant. 

For t > the time evolution is governed by the Langevin equation 

where r/(x,t) is the gaussian white noise produced by the thermal bath 
at the temperature of the quench T, jj, is the set of parameters entering 



the free energy functional 7i[$, p\ and p = for NCOP while p = 2 for 
COP. 

In the following we take a free energy functional of the Ginzburg- 
Landau form 

d = / d d x [|(V$) 2 + ^$ 2 + |($ 2 ) 2 ] (6) 

with /i = (r, g, i?). The equation of motion then becomes 

= -T(-iV) p [-i?V 2 $ + r$ + g&] + 77(f, *). (7) 

In the study of deep temperature quenches usually one puts T = (i.e. 
T)(x,t) = in eq. (0)) regarding the temperature as an irrelevant param- 
eter and (r, g) are chosen in the sector r < 0, g > corresponding to 
final equilibrium states inside the ordering region. However, as it will 
be clear shortly, it is of considerable interest also the edge of this sector 
(i.e. r = 0, g > 0), even though no phase ordering occurs there since 
the equilibrium value of the order parameter vanishes. Thus, the states 
relevant to our discussion are those in the T = plane of the (T, r, g) 
space with r < and g > 0. We will set R — 1 and T = 1 for simplicity. 

3 Numerical results 

As stated in the Introduction, the quantity of interest is the equal time 
order parameter correlation function 

G(\x-x'\,t) =< $(x,t)$(f ,t) > - < $(x,t) > 2 . (8) 



We consider first the quench to the trivial state with //j = (r = 
0, g = 0). In this case eq. (0) can be solved exactly and one finds 

C(k,t)=Ae- 2k2+n . (9) 

Therefore, defining L(t) = (2t) 1/{2+p) the real space correlation function 
is in the scaling form (1) with a — d, z — 2 for NCOP, z = 4 for COP 
and F(x) = exp(— x 2+p ). Notice that in this case eq.(9) is not just 
an asymptotic behaviour, but it is obeyed exactly along the whole time 
history, from the instant of the quench onward. In the language of critical 
phenomena this means that the width of the critical region is maximally 
amplified, warranting the identification of \L\ with a (trivial) fixed point 
on the T = manifold. The next step is the exploration of the domain 
of attraction of this fixed point, the search for other fixed points and for 
the crossover induced by their competition. 

Since for g ^ the theory is not soluble we proceed by numerical sim- 
ulation. If scaling holds, from eq.(l) we have S{t) = G(0,t) ~ L~ a (t) ~ 
t~ a / z . In the following we will use the behaviours of L(t) and S(t) to get 
the pair of exponents (a,z). The evolution of G(\x — x'\,t) is obtained 
by numerical solution of eq. ([7]), for a two-dimensional 100x100 lattice. 
The equal time correlation function is obtained by averaging over dif- 
ferent realizations of the time histories (the number of such realizations 
ranges from 1 to 20, according to the quality of the data). We have found 
convenient to extract the characteristic length L(t) from the half-height 
width of G(\x — x'\, t). 
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We consider first the behaviour of eq. (0) for quenches on the r = 
axis. In this case no phase ordering occurs, in the usual sense, since even- 
tually the order parameter vanishes. Nevertheless one can still observe 
the formation and subsequent growth of domains which can be defined by 
locating their boundary on the contour £) = 0. L(t) can be thought 
as the typical length associated to these structures. 

The results of our simulations are presented in fig. 2 for NCOP, where 
the behaviours of L(t) and of S(t) are displayed for different values of g. 
We do not observe significant differences among the behaviours of L(t) 
as g is varied and we conclude that Lit) obeys an asymptotic growth 
law with an exponent consistent with z = 2. The quantity Sit), on the 
other hand, shows a more complex behaviour, in that the asymptotic 
scaling S{t) ~ t~ a l z sets in after an initial transient, which widens as g 
gets bigger. For g — 10 the value of the exponent a/z — 1.2 is slightly 
different from the one found for g — 0.1 and g = 1, where a/z = 1 
implying a = 2. This is due the initial transient which for g = 10 is not 
completed over the time of the simulation. 

In conclusion the data show an asymptotic scaling behaviour identical 
to the one found in the quench to the trivial fixed point with z = 2 and 
a = d. The critical regime where scaling holds shrinks moving away from 
the trivial fixed point along the r = axis. 

The analogous results for COP are presented in fig. 3. Again scaling 
behaviour is obeyed by L{t) and S{t) with exponents consistent with 



8 



z = 4 and a = 2 but, differently from NCOP, this feature sets in almost 
immediately after the quench, for every value of g. Both for NCOP and 
COP quenches to the r = axis are controlled by the trivial fixed point. 
Therefore, all quenches with r = fall into the same universality class 
and the cubic term in eq. (0) is asymptotically irrelevant. 

We turn, now, to the phase ordering region by setting r < and 
g > 0. When considering quenches inside this region, as it has been 
well documented in the literature [|J, the asymptotic scaling behaviour 
is different from the one just found on the r = axis. In the phase 
ordering region there exists a non trivial fixed point characterized by 
(z = 2, a = 0) for NCOP and (z = 3, a = 0) for COP, whose domain of 
attraction is the entire sector (r < 0, g > 0). For quenches sufficiently 
close to r = axis, but inside the ordering region, we expect to observe 
the crossover from trivial to non trivial scaling behaviours. 

We have performed simulations with g = 1 and three different values 
of r, both for NCOP (fig.4) and COP (fig.5). For NCOP, L(t) scales with 
the same exponent z = 2 both when the dynamics is dominated by the 
fixed point of the ordered region or by the trivial one. More interesting 
is the behaviour of S(t), which is characterized by the sequence of three 
regimes. In the early regime, whose duration is shorter the farther away 
from the r = axis the quench occurs, we find that S(t) decreases with 
a behaviour similar to the one found in the case of a quench on the 
r = axis. During this early regime the order parameter evolves as 
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if the potential did not have a double well and relaxes locally toward 
zero in order to eliminate spatial inhomogeneities as in the quenches to 
the trivial fixed point. This occurs because, as a consequence of the 
initial disordered state, the gradient term dominates in eq. (y). Then, at 
some time t\ (see the schematic representation of fig.l), the double well 
structure starts to play a role. At this point S(t) stops decreasing and 
enters the intermediate time regime characterized by exponential growth 
of the order parameter toward local equilibrium at the bottom of the 
wells. In this regime L(t) is approximately constant, as in the linear 
theories ||. The intermediate regime terminates with the formation of 
domains within which the order parameter is close to saturation. From 
this point onward the late stage is entered with dynamics dominated by 
interface motion. Since the system now is close to saturation S(t) ~ 
S'(oo) = — |, implying scaling controlled by the non trivial fixed point 
with a = 0. 

The case of COP, shown in fig. 5, is particularly interesting because 
this whole crossover structure is manifested also in the growth law of 
the size of domains L(t). In the first regime, whose duration is again 
controlled by the distance from the r = axis, the dynamics is dominated 
by the trivial fixed point, with scaling exponents z = 4 and a = d. As 
stressed before for quenches on the r = axis, no initial transient is 
observed for COP. This allows a precise determination of the exponents 
a and z in the early regime obtaining a = 2 and z = 4. Then, when the 
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system feels the presence of the local potential, the exponential growth of 
S(t) is accompanied by an approximately constant behaviour of L(t) (or, 
at least, by a lower growth rate) which characterizes the intermediate 
time regime. Later on, when S(t) saturates to its asymptotic value, the 
system enters the late stage, dominated by the non-trivial fixed point 
with a = and z = 3 (this is not shown in fig. 5 since the z exponent 
reaches its asymptotic value for much longer times). 

Finally we consider the effect of performing off-critical quenches, in 
the case of a conserved order parameter, by varying m =< $(x, t) > 
through the region m 2 < —- g (see fig. 6). The asymptotic scaling be- 
haviour inside this region is still controlled || by the same non trivial 
fixed point of the symmetric quenches, with z = 3 and a = 0. The 
domain of attraction of this fixed point is, therefore, the entire region 
m 2 < — - where phase separation occurs. At the intersection of the co- 
existence curve with the T = manifold (i.e. at the very edge of this 
region, with m 2 = — £) no phase separation occurs since the field reaches 
a uniform configuration. In this case the amplitude of G(\x — x'\,t), 
namely S(t), vanishes asymptotically approaching equilibrium, as for a 
critical quench on the r = axis. More precisely in this point of the 
phase diagram the long time behaviour of the model is no longer con- 
trolled by the same fixed point as for critical quenches, which on the 
contrary requires S(t) — > constant, and we expect to observe a dynamics 
reminiscent of the r = 0, m = case, where also S(t) — > 0. Therefore 
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it is of considerable interest to understand the nature of the quenches 
close to m 2 = — ^ and to study the crossover phenomena induced on the 
dynamics of the system for m 2 < ~. 

In order to address these questions we have simulated asymmetric 
quenches by preparing the initial configurations with different values of 
m 2 in the range (0,— p. First of all it must be stressed that, for quenches 
inside the metastability region, i.e. for m 2 sufficiently close but differ- 
ent from — ^ (see fig. 6), the influence of the ordering fixed point with 
(a = and z = 3) is hardly observed, because the system relaxes into 
a metastable state inside a single minimum of the local potential (the 
one of the majority phase). This can be avoided by increasing the vari- 
ance A of the initial condition. When stable equilibrium is reached in 
the final state, the behaviour of S(t) reveals again the presence of three 
distinct dynamical regimes. In fig. 7 this quantity is shown for quenches 
to (r = — 1, g — 1) with different values of m. In the early regime, 
whose duration increases indefinitely approaching the coexistence curve, 
we again find the trivial scaling behaviour with (a — d, z — 4) as for 
r = and m — 0. Later on the system enters the intermediate and 
eventually the asymptotic regime, with the same features described in 
the case of critical quenches. In the renormalization group language this 
means that, by increasing m from toward the coexistence line, the tra- 
jectories start closer to the domain of attraction of the trivial fixed point 
r = 0, m = 0. This is analogous to the amplification of the early regime 
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in the critical quench by decreasing the value of r. However, asymptoti- 
cally, the non-trivial fixed point prevails for m 2 < — ^. Therefore in the 
space of the parameters (r, g, m) the region m 2 = — | and the axis r = 
axis are the domain of attraction of the trivial fixed point. 

We conclude this section observing that, as far as the role of the 
parameter R in eq. (0) is concerned, by rescaling one can easily show 
that the effect of increasing R corresponds to a magnification of the 
space and time scales. Since R controls the range of the interaction it is 
conceivable that the early stage scaling regime is more clearly observable 
in systems with sufficiently long range interactions. 

4 Comparison with the large- N model 

In this section we compare the results of our simulations with the solution 
of the large- N model . When the order parameter is an ^-component 
vector field = (<&i(x), $jv(x)) in the limit of an infinite number 
of components eq.(7) is linearized 



From the exact solution of this model one finds that the asymptotic scal- 
ing properties depend as in the scalar case on the pair of coupling con- 
stants r and g. More specifically one finds that there is a universality 
class, under each heading NCOP or COP, for each of the following three 
regions in the (r, g) space 




(10) 



where S(t) must be computed self-consistently through S(t) =< $ 2 



(x,t) >. 
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/ii = (r = 0, g — 0) (trivial critical state) 

/jl 2 = (r — 0, g > 0) (non trivial critical states) 

/i 3 = (r < 0, g > 0) (phase ordering region). 

In the renormalization group language this means that there are three 
fixed points and the extension of the universality classes depends on the 
relative stability of these fixed points. Actually, for quenches to fi 2 there 
is a critical dimensionality d c , which depends on the initial condition, 
above which the non linearity in the problem becomes irrelevant. For an 
initial condition of the type 

c(M = o) = p (ii) 

one finds d c = di + 9, where di is the lower critical dimensionality of the 
static problem. With 6 = one has d c = di = 2. Therefore, solving the 
model for d > 2 one finds the following asymptotic behaviours for NCOP 

C(k,t)~F(x) (12) 

/i 2 ) C(k,t) ~ F(x) (13) 

/i 3 ) C(M) ~ L d (t)F(x) (14) 

with L(t) ~ t 1 ' 2 and F(x) = e'* 2 . 

Let us briefly comment on these results. First of all we have z = 2 
everywhere. For quenches to the trivial state \i\ the large- iV model and 
the scalar model coincide since the local potential is absent. The same 
expression (9) and (12) is found for the structure factor. For quenches 
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to /j,2 the long time behaviour is the same as for the quench to /ii, up to 
corrections to scaling, as expected since d > d c . The scaling behaviour 
is characterized by (a — d, z — 2) both for quenches to H\ and /x 2 . 

Different is the case of quenches inside the region of coexisting phases, 



where eq. (|14| ) is asymptotically obeyed in any dimension, showing that 
no upper critical dimensionality exists and therefore that the non linear- 
ity of the problem is always relevant. In this case scaling is characterized 
by (a = 0, z = 2). Furthermore, computing the early time behaviour 
for \r\ sufficiently small a scaling regime is found which is identical to 
the one found for quenches to This behaviour is separated from the 
asymptotic scaling regime described above by an intermediate regime of 
exponential growth. 

In the large-iV model the variety of asymptotic properties is more 
complex when quenches with COP are considered, due to the existence 
of multiscaling. The following behaviours are found 

/n) C(k,t)~F(x) (15) 

H 2 ) C(k,t) ~ F(x) (16) 

Us) C(k,t) ~ L(t) a M (17) 

where L(t) ~ t 1 / 4 , F(x) = e~^ xi+c ^ and a(x) = d[l — (x 2 — l) 2 ] is a 
generalized x-dependent exponent. The general discussion of fixed points 
and their relative stability goes along the same lines as for NCOP. 

Now we make a qualitative comparison with the results for the scalar 
case presented in the previous section. The first observation concerns 
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the existence of an early stage scaling regime which is observed in both 
models for quenches in the phase ordering region, with the same expo- 
nents a and z. This was to be expected since the very nature of this 
behaviour is due to the trivial fixed point with r = g = 0, where the 
two models coincide. We emphasize that the occurrence of early scaling 
is particularly interesting for conserved scalar order parameter since in 
this case, differently from the large- N limit, a crossover from z = 4 to 
z = 3 exists in the growth law. Secondly, we comment on the existence 
of a critical dimensionality for quenches to ji2- In the large- N model 
d c = di = 2. In the scalar case d\ = 1. However, we don't know whether 
d c is still equal to d t . The results of section 3 indicate that d = 2 is above 
the critical dimensionality. For what concerns the late stage regime, the 
two models are quite similar for NCOP, since scaling is obeyed with the 
same exponents. With COP, on the contrary, there is less of an analogy 
since the vector model exhibits multiscaling. 

5 Conclusions 

In this paper we have analyzed the quench to zero temperature of a sys- 
tem with a scalar order parameter through the numerical solution of the 
time dependent Ginzburg-Landau equation, both with conserved and 
non conserved order parameter. We have paid particular attention to 
the early stage of the quench. According to the values of the parameters 
(r, g) which characterize the final equilibrium state, an early stage scaling 
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regime associated to pure diffusive behaviour may be observed before the 
usual late stage scaling regime sets in. This is a crossover pattern of the 
same type found in the analytical solution of the large- iV model. In that 
case the pure diffusive behaviour is associated to the trivial fixed point 
(r = 0, g — 0) while the late stage scaling behaviour is associated with 
the non trivial fixed point (r < 0, g > 0) of phase ordering. The competi- 
tion between the two fixed points determines the crossover pattern. The 
question of the observability of this phenomenon is clearly related to the 
time span of the early stage. In the paper we have mentioned long range 
forces and off critical quenches toward the coexistence curve as means of 
amplification of the early stage. 
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6 Figure captions 

Fig.l - Schematic representation of the dynamical response of the scalar 
model, both for NCOP and COP. The three regimes are shown, separated 
by the crossover times t\ and t 2 . 

Fig. 2 - Behaviour of L(t) (fig. 2a) and of S(t) (fig.2b) for a non con- 
served system quenched on the r = axis, with different values of g 
(g = 0.1, g — 1 and g = 10). The continuous lines represent respectively 
the power laws t\ and t~ x associated to trivial scaling. Best fit yields 
z = 2.01 ± 0.05 for every value of g and a/z — 1.0 ± 0.1 for g — 0.1 and 
= 1, a/z = 1.2 ±0.1 for g = 10. 

Fig. 3 - Behaviour of L(t) (fig. 3a) and of S(t) (fig. 3b) for a conserved 
system quenched on the r = axis, with different values of g (g = 0.1, 
g = 1 and g = 10). The continuous lines represent respectively the power 
laws U and t~z associated to trivial scaling. Best fit yields z = 3.98±0.05 
and a/z = 0.51 ± 0.03, for every value of g. 

Fig. 4 - Behaviour of S(t) for a non conserved system quenched inside 
the phase consistence region, with g = 1 and different values of r (r = 
—0.005, r = —0.05 and r = — 1). The duration of the early scaling regime 
increases as \r\ decreases. 

Fig. 5 - Behaviour of L(t) (fig. 5a) and of S(t) (fig. 5b) for a conserved 
system quenched inside the phase consistence region, with g = 1 and 
different values of r (r = —0.03, r = —0.3 and r = —1.3). The continuous 
lines represent the trivial scaling behaviour with z = 4 and a = 2. Best 
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fit of data with r = -0.03 yields z = 4.1 ± 0.1 and a/z = 0.49 ± 0.05. 

Fig. 6 - The phase diagram in the (m,T) plane. The phase coexis- 
tence region and the spinodal line, above which metastability occurs, are 
shown. 

Fig. 7 - Behaviour of S(t) for a conserved system with r = —1, g — 1 
and m ranging between and 0.7 (from the top to the bottom: m = 0, 
m = 0.1, m = 0.2, m = 0.3, m = 0.4, m = 0.5, m = 0.6 and m = 0.7). 
The lower straight line represents trivial scaling with a/z = 0.5. 
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